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In earlier papers [11-17], it is found that the Ricci scalar behaves in a dual 
manner (i) like a matter field and (ii) like a geometrical field. In ref. [17], in- 
homogeneous cosmological models are derived using dual roles of the Ricci 
scalar. The essential features of these models is capability of these to ex- 
hibit gravitational effect of compact objects also in an expanding universe. 
Here, production of spinless and spin- 1/2 particles are demonstrated in these 
models. 
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1. Introduction 



Consideration of higher-derivative gravity is not new and it is a strong 
candidate in the context of string theories and Einstein-Gauss-Bonnet grav- 
ity. In the context of the early universe, use of higher-derivative gravity 
started from more than three decades ago. In 1977, Stelle considered it and 
pointed out that, in contrast to square of Weyl tensor, introduction of R 2 
or f(R) terms did not lead to the ghost problem [1]. Here R is the Ricci 
scalar and f(R) is the linear combination of R 2 terms subject to the condi- 
tion lim R ^ f(R)/R = [2]. From these theories, it is found that the Ricci 
scalar behaves as a physical field also [3-10]. 

While quantizing gravity (quantizing components of the metric tensor), 
this theory has problem at the perturbation level where ghost terms appear 
in the Feynman propagator of gravitons. Around 25 years back, Starobinsky 
noted that if coupling constants in the action of higher-derivative gravity 
are taken properly, ghost terms do not appear and only one scalar particle, 
represented by the Ricci scalar, is obtained with positive energy as well as 
positive squared mass [6,7]. He called this particle as "seal aron" [6,7]. In pa- 
pers [1-8], gravitational constant G is either taken equal to unity or lagrangian 
density is taken as j^q(R+ higher - derivative terms). As a result, (mass) 2 
of the Ricci scalar, does not depend on the gravitational constant G. As 
G has a very important role in gravity, in papers [11 - 17] as well as here, 
lagrangian density is taken as ( 1 q] tG R+ higher - derivative terms) leading to 
a drastic change where (mass) 2 of R depends on G also. In these papers [11 - 
17], choosing coupling constants of the higher-derivative gravitational action 
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in a suitable manner, it is found that the Ricci scalar behaves like a physical 
particle, called as " riccion", with positive (mass) 2 without any ghost prob- 
lem. This feature leads to manifestation of Ricci scalar in dual manner (like a 
geometrical field as well as a scalar particle). It is important to mention that 
riccion is a scalar particle (having only one mode unlike a graviton having 
five modes), but it is different from scalar mode of graviton. Detailed discus- 
sion on this difference is given in Appendix A. Riccion is given by a scalar 
field R = r]R, where rj is a constant having length dimension. A "riccion" is 
different from a "scalaron" in two ways (i) mass dimension of riccion is one 
like other scalar fields used in a physical theory, such as quintessence, infla- 
ton and Higg's particle, whereas mass dimension of scalaron is two and (ii) 
(mass) 2 of scalaron does not depend on G, whereas riccion (mass) 2 depends 
on G. 

In refs. [11-17], some interesting results are obtained using dual role of the 
Ricci scalar. One of the important consequences of the same is derivation 
of models of the early universe using a physical method of phase transition 
and spontaneous symmetry breaking. In ref. [11], spatially homogeneous cos- 
mological models were derived using these methods and it was shown that 
universe should bounce at the critical temperature ~ 10 18 GeV. In ref. [15], 
employing these methods, models of the early universe are derived exhibiting 
inhomogeneity and anisotropy at small scales. It is found that these models 
are asymptotically homogeneous and isotropic. 

The theory of gravity is manifested through geometry of the space-time. 
Two types of space-times are used for manifestation of gravitational interac- 
tions. One type is static space-times representing geometry around compact 
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objects e.g., Schwarzschildspace-time, Riessner- Nordstrom space-time and 
others. The other type of space-times are homogeneous and inhomogeneous 
models representing the expanding universe. Former type of models do not 
give any information about the expanding universe and the latter type of 
models are unable to account for gravitational effect of individual objects 
in the expanding universe. So, models exhibiting local gravitational effect 
of an individual object as well as expansion of the universe simultaneously 
are required to discuss evolution of the universe in a more natural way. Us- 
ing dual roles of the Ricci scalar, such cosmological models are derived in 
the reference [17], giving the first two stages of the beginning of the early 
universe. 

In the present paper, production of spinless and spin-1/2 particles is 
demonstrated in the inhomogeneous models derived in the reference [17] . 
Contributions of these particles to the development of the universe will be 
discussed in paper II. 

The paper is organized as follows. For convenience of the reader, a review 
of reference[17] is given in section 2. 

In section 3, taking interaction of riccions with another scalar field <3>, 
production of spinless particles is demonstrated. In section 4, creation of 
spin-1/2 particles has been shown. Cosmological consequences of these cre- 
ated particles are planned to be discussed in the second paper of the series. 

Natural units (h — c — 1) are used throughout the paper. Here H and 
c have their usual meaning. For the sake of convenience, it is given that 
lGeV = 1.16 x 10 13 K = 1.78 x 1(T 24 g, 1 GeV^ 1 = 1.97 x 1(T 14 cm = 
6.58 x 10" 25 Sec. 
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2. (a) Dual role of the Ricci scalar and riccion 



In this section, a brief review of the reference [17] is given. Here, the 
gravitational action with thermal correction is taken as 



where a and A are dimensionless coupling constants. Here T is the temper- 
ature, G is the gravitational constant and rj is another constant with length 
dimension, which is used for dimensional corrections only [11-17] . The in- 
variance of S g , given by eq.(2.1) , under transformations g^ v — > g^ v + 5g^ u 
yields the gravitational field equations 




R 



16r] 2 \a{R 3 + 9UR 2 ) + ^-R 3 ] , 



(2.1) 



■ ~ 16^G + 4AaT2 l ^ ~ 2 9ilvK) + a( - 2R ^ ~ 2 9^ DR 
-\g^R 2 + 2RR»u) - 16 V 2 \a{3Rl u - 3g^UR 2 

+H-\g^OR 2 + 2R^UR + R 2 ^ - X - gilv R? + 3R 2 R ilv )} = 0, (2.2) 



where 




(2.3) 



General anastaz for space-time is given as 



dS 2 = g^dx^dx". 



(2.4) 



5 



Trace of field equations (2.3) yields 
1 



1 



^ + 2XaT 2 R + a(8DR + -R 2 ) + 8r] 2 \aR 3 = 0. (2.5) 
327rG j 2 



i?, being a linear combination of second-order derivatives and square of 
first-order derivatives of g^ u , has mass dimension 2, whereas the same for 
scalar fields in known theories is 1. So, eq.(2.5) is multiplied by i] (having 
length dimension) and r]R is recognized R. Thus R, having mass dimension 
1, represents a scalar particle. This scalar particle is called "riccion". In 
the case of scalaron, r\ is dimensionless and its magnitude is 1 [6]. So, it is 
obtained that 



1 



+ 2XT' 2 R + Xi] 2 R 3 = 



(2.6) 



32irGa 

If R is a basic physical field, there should be an action yielding eq.(2.6) 
for invariance of S 1 ^., under transformations R — > R + 5R. 

In what follows, is obtained. If such an action exists, one can write 



5S 



R 



- d A xV^SR[(a + —RR+- 



32nGa 



+2AT Z 



R+Xr] 2 R s ] (2.7a) 
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which yields eq.(2.6) if SS R = under transformations R — > R + 5R. 
Eq.(2.7a) is re- written as 



5S 



r = J d 4 x^j5 



-d^RdJSR) - —RR 2 - V T (R) 
2 MV ; 48?7 v ; 



(2.76) 



where 



V T (R) 



l -m 2 (R 2 + i-T 2 ) + \XR A + \XT 2 R 2 - %-T A 



l 

12 



1 

I' 



90' 



(2.7c) 



with 
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R, R and d 4 Xy^g are invariant under co-ordinate transformations. So, 
R(x) = R(X),R(x) = R(X) and d 4 x^g = d 4 X, where X\i = 0,1,2,3) 
are local and x l {i = 0, 1, 2, 3) are global coordinates. Moreover, 
d 2 1 _do mri .... d dq ij d d 2 



a = g 



_ , 1 mn d 9mn jj dg^_8_ = 

dx*dx3 2 y dx i y dxi dx i dxi dX*dX3 



in a locally inertial co-ordinate system, where gij = rjij (components of 
Minkowskian metric) and g L \ = (comma (,) stands for partial derivative). 
Thus, in a locally inertial co-ordinate system, 



5S 



R 



d 4 X5 



-'I 



d 4 X 



\d»Rd»{5R) - -±-RR 2 - V T (R) 
l -d»Rd,{5R) - -^RR 2 - V T {R) 



Employing principles of covariance and equivalence in eq.(2.8) 



5S 



k = 5j d 4 x{. 



-d^RdJSR) - —RR 2 - V T (R) 



yielding the action for riccion as 



(2.8) 



S 



R 



d 4 x^—g 



1 



-d»Rd,(5R) - -RR 2 - V T (R) 



(2.9) 



It is important to mention here that R is different from other scalar fields 
due to dependence of (mass) 2 for R on the gravitational constant and the 
coupling constant a which is given by the eq.(2.1). Moreover, it emerges 
from geometry of the space-time. 



Further, action for a scalar and a Dirac spinor ip are taken as 

S(m) = J d A X< 

+ ^ J d 4 xy^[i>i^D^ - aR^ + c.c], 



(2.10) 



where A and a are dimensionless coupling constants. Here $ = ^7°, are 
covariant derivatives. 7 M are curved space Dirac matrices and 7^ being flat 
space Dirac matrices. 
So, the total action is 



S R + S rn 



d^Rd^SR) -^RR 2 - V T {R) 



+ J d*x{V=g[\<rd li $d v & - Ur 2 ^ 



+ 7; / d^Xy^glipiYD^ip - crRipip + c.c.}, 



(b) Models of the early universe inspired by dual role 
of the Ricci scalar 



(2.11) 



Here models of the early universe are obtained in vacuum states of riccion, 
given by Higgs like potential (2.7c) [11, 17, 18] employing the condition 



T 



dV 
~dR 



0, 



8 



or 




(2.12) 



This equation yields turning points of V T (R) as 



R = 



(2.13a) 



and 



R = ±-/(T 2 -T 2 ), 



(2.136) 



where 




(2.13c) 



These states can be written in terms of vacuum expectation value of R as 



the field R remains confined to the state < R >= 0. But when temperature 
goes down such that T <^ T c , R tunnels through the temperature barrier 
T = T c and settles in the state < R >= (1/2)T C when T <C T c . The state 
< R >= is called the false vacuum and < R >= (1/2)T C corresponds to the 
true vacuum state . Moreover , symmetry is intact in the state < R >= 0, 
but it is broken spontaneously at < R >= (1/2)T C and energy is released as 
latent heat in the form of radiation with density given as 



The line - element for the cosmological model of the early universe is 
taken as 



< R >= and < R >= (l/2)^/(T c 2 — T 2 ). It shows that as long as T > T c 




m 



4 



4A 



(2.14) 



dS 2 = dt 2 - a 2 {t)f{r) dr 2 + r 2 {dd 2 + sm 2 ## 2 ) , 



(2.15a) 
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where t is the cosmic time ,0 < r < oo, < 6 < n and < ip < 2tt. 

Looking at the line - element , given by (2.15a) , one may think that the 
function / can be absorbed by redefining r as new coordinates f such that 

f(r) [dr 2 + r 2 {d6 2 + sin 2 6d^ 2 )} = dr 2 + r 2 {d® 2 + sin 2 Qd^ 2 ) (2.156). 

It is important to mention here that redefinition of coordinates are not valid 
at points where / = or oo, as such points lead to spatial singularity 
(singularity at points on t = constant hypersurface) . So,the absorption of 
/ in r coordinates is possible only when / ^ and / < oo at all points of 
t = constant hypersurface. It means that only after getting explicit definition 
of / , it is possible to decide about /. This aspect will be looked into below. 

Using geometrical aspect of the Ricci scalar , one can compute R in the 
background geometry of the line - element , given by eq.(2.15) , as 

p = \«( a i n 2 i 1 \ 2 ( d2f2 1 2g/2 i i 3 ( df2 



dr 



2i 



(2.16) 



a*r* L f 4 ^ dr 2 r dr > 2/ e 

which is linear sum of the homogeneous part , given as 6??j^ + (^j | as 
well as the inhomogeneous part containing derivatives of f 2 . < R > , being 
vacuum expectation value of R , is homogeneous . So , in the vacuum states 
< R >= and < R >= ±(1/2)T C , 



a?r 2 



2 ( d 2 f 2 2df 



-{ 



+ 



3 ( df 



First stage of the early universe 



} + ^( 



2/6 v d r 



0, 



(2.17) 



Coming to the physical aspect of R, one finds that at T > T c , R stays in 
the false vacuum state R =< R >= 0. So , from eq.(2.16) , one obtains the 
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differential equation 

which yields the solution [17] 

a 2 = a 2 c + -L^, (2.19) 

where t p i is the Planck time. At this point, a question arises whether a c = 
or it is non - zero . In what follows, this question is answered . 

Components of energy - momentum tensor for matter field R can be 
defined as [11, 14-15] 

rp Q 1 SS R 

which is obtained as 



T„ v = d^Rd u R - g^l-ffiRd^R - V T (R)] - 2r]R^DR + 2r ] m 2 RR l 

-2r)\R 3 R^ - -rj\T 2 RR^ + 4 V m 2 R.^ - 4 V m 2 g^nR 

-47 7 A J Rf^ + AriXg^DR - \T 2 R^ + \T 2 g^UR. (2.20) 
In the false vacuum state R =< R >= 0, one obtains [11,16, 17] 

T^UW < 0, (2.21) 

where are given by eq.(2.20) and are velocities normalized as 

= +1. 

The result, given by the inequality (2. 21), shows that the energy condition is 
voilated in the false vacuum state when T > T c [11,15,17]. It means that 
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the universe will bounce without any encounter with temporal singularity 
(singularity in time). As a result , 

a c ^ 0. (2.22) 
The partial differential equation (2.17) yields an exact solution 

f= h _I°1 4/7 . (2.23) 

L r J 

Thus, in the false vacuum state, model of the early universe is obtained 
as given by eq.(2.15a) with a 2 (t) and f 2 (r) defined by eqs.(2.19) and (2.23) 
respectively. 

The expansion of the cosmological model, derived above, is adiabatic. 
So, entropy will remain conserved , which implies that the average uniform 
temperature will fall as 

r n-1/2 

T(x\a 2 c + (| t \/t p l)\ . (2.24) 
Second stage of the early universe 



When T < T c 

R=<R >= ±(1/2)T C . 
In this state also , one obtains two differential equations 

<+©V4 T < (2 - 25) 

and the other differential equation is the equation given by eq.(2.17). Eq.(2.25) 
yields two cases. The first case is 



rd /a\2i 1 
6 "{a + U I = 2 T " (2 - 26) 
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yielding the solution 



a 2 (t) 



I 24 



sink 



(t-t u )J^ + D 



a 2 e [sinhD] 1 sinh (t — t le )\j -+■ D 



a le exp 



(2.27a) 



with ai e = a(t = ti e ) and 



D = sink 1 



"lei 



24 



Here ti e is the time at the end of first stage of the universe and a± e is the 
corresponding scale factor, given as. 



a(t) ~ ai e exp (t — t le 



2Ar] 



The second case is 



6 "{- + (-) 2 } = -^ 



(2.276) 



(2.28) 



which integrates to 



r n 
a 2 (t) = a\ Q sin\{t - t )J ^ + - 



(2.29) 



In true vacuum states also ,/ is given by eq.(2.23). 

But a(t), given by eqs.(2.27), shows that the cosmological model will 
expand more rapidly in the state R =< R >= (1/2)T C compared to the state 
R =< R >= 0. The expansion , in the state R =< R >= — (1/2)T C , is 



13 



given by eq.(2.29), which shows that — a\ < a 2 < a\ Q leading to imaginary 
a(t) also, which is unphysical. So the state < R >= — (1/2)T C will not be 
considered any more. 

Thus one finds that vacuum states R =< R >= and R =< R >= 
(1/2)T C lead to cosmological models 



dS 2 = di 2 - 



tpi 



1 - 



r 



4/7 



dr 2 + r 2 (dd 2 + sin 2 9d^j 2 ) (2.30) 



at temperature T > T c and 



dS 2 = dt 2 — a\ e exp 



(*-*le)l 



ro 
r 



4/7 



dr 2 + r 2 {d6 2 + sin 2 9di) 2 ) (2.31) 



when T < T c . 

Divergence of scalar polynomials of the curvature tensors at points r = 
and r = r . shows that models ( given by eqs.(2.30) and (2.31)) are singular 
at these points [21-22]. Moreover, these models are inhomogeneous at small 
scales i.e. when r is comparable to r and homogeneous when r >> r . 

The model (2.31) exhibits accelerated expansion of the universe which 
is consistent with consequences of recent experiments la Supernova [24] and 
WMAP [20]. 

It is interesting to note that models, given by eqs.(2.30) and (2.31), reduce 
to homogeneous Robertson- Walker type models as r — > oo. As far as spatial 
singularity is concerned, at a particular instant of time , these models can be 
compared with Scwarzschild space-time [19, 21] 



dS 2 



1 - 



2GM 



dt 2 - 



1 - 



2GM 



dr — r z 



d6 2 + sinHdnf 



(2.34) 



It is remarked that the Scwarzschild space-time given by eq.(2.34)is singu- 
lar at r = and r = 2GM, but only r = is the real physical singularity and 



14 



the other point is a co- ordinate singularity. In models , given by eqs.(2.30) 
and (2.31), both points r = and r = r exhibit real singularities as men- 
tioned above. Physically r signifies radius of particles or compact objects 
dominating the expanding universe. Further investigations are carried out 
for r > r . 



3. Creation of spinless particles 

Invariance of S, given by eq.(2.11) under transformation <3> — > $ + 5$ 
leads to the Klein - Gordon equation 

[□ + M 2 ]$ = 0, (3.1) 

where 

M 2 = A < R > 2 . (3.2) 

In the background geometry of the line - element (2.15 a) with eq.(2.23), 
the equation (3.1) is re-written as 



<9 2 $ 3a <9$ 
~W + ^~8t 
1 



1 - 



r 



-4/7 d 2 & 
dr 2 



r 2 a 2 



1 - 



r 1 -6/7 8 
dr 



9$ 

dr 



a 2 r 2 





-4/7 


L r _ 





1 d / . . 9 x 1 



<9 2 



sm 2 ^ <9-?/> 2 



$ + M 2 $ = 0, (3.3) 



-sin9 09^ 89' 

where dot (.) over the variable denotes the derivative with respect to time. 
The solution for this equation is taken as 



(27r)VV(l - -f) ]" £ [^ m ^ m (0^ m (^V)e^ + c.c], (3.4) 

where r > r , m = —/••• + /,/ = 1, 2, 3, • • • and — oo < A; < oo. Yi m (9, ip) 
satisfies the equation 

i d 2 



' ' sin9^-\ + 



sin9 89 



89' sin 2 9 dip 2 
15 



Yim — —1(1 + l)Yi m 



(3.5a) 



with normalization as 



J sin9d9dipYi m Y llm , = S w S r , 



(3.56) 



Connecting eq.(3.3) with eq.(3.4) and using eq.(3.5), it is obtained that 

-4/7 



Ak.r 



A + 3 % I (( k2 | *(* + *) | 6r (2r-r 



a 2 r 2 7a 2 r 2 (r — tq) 2 



= 0. 



(3.6) 



Using the convolution theorem [25] and integrating over r, eq.(3.6) looks 



like 



[°° e lk - r dr\^ klm + —^ klm + {4- / 



k 2 r°° 



-ik.y 



1(1+1) 



-ik.y 



a 2 T] Jro+e y 2 



1-^ 



-4/7 



dy + 



T] Jro+e 
6r roo 



1 - 



r 



-4/7 



(2y - r ) 
7a 2 r)Jr +ey 2 (y - r ) 2 



1-1° 



-4/7 



+Ml}m klm 



o. 



dy 
(3.7) 



Here e is an extremely small positive real number. 

Details of the evaluation of integrals with respect to y are given in the 
Appendix B. Using these results, in eq.(3.7), one obtains 

rX k 



T 3fl • 

a 



^klm + M 2 



fcim 



o, 



(3.8a) 



where 



X klm = (l/a 2 r ? )e 3 / 7 (r + ef' 7 cos{k(r Q + e)} 



fc2+ /(/ + l) 



(r + e) 



-yr e- 2 (r + e)" 1 + ^r V 2 (r + e)" 



(3.86) 



Case 1: The case of state < R >= 
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Using a(t) from the line-element (2.30) for this state, the equation (3.8) 
is written as 



klm^P 



^klm = 0, 



which yields the solution for t > as 

Vkim = r 1/2 [AJ-i/ 2 (r) + BY_ 1/2 (r) 



where A and B are integration constants. Here r is defined as 

T = \JtpXkim(t + tpal ) 

showing that r — > oo when £ — > oo. 
For large r, 

cos(r) sot(t) 
J-i/ 2 (t) ~ and F_ 1/2 (r) ~ 



vrr/2 



vrr/2 



So, when r is large, 



#fczm = [7rr 2 /2] l/2 [Acosr + Bsinr] 

' r [l + i)e~ %T + (I - iY 



i 



2t P X k i m 



1/2 

r 



using the normalization condition 



{^klm^klm) = 1 = -(^fc, ro ,^fc, m ), 



where 



®klm — 



(2nrj- 1 W) 1 / 2 r(l- T f) 1/7 
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^kim{t)e lk - r Y lm (3.126) 



(3.9) 



(3.10a) 



(3.106) 



(3.11) 



(3.12a) 



and the scalar product is defined as 

"CO /"7T p2lT 



{®kim,®kim) = ~i I I I V=g^dZ[<$> k i m (dt<f>* k/l , m ,) 

Jr +eJ0 JO 



-{d t ® klm )®l, llm ]\Y lrn {6M 2 , 

where E is the t — t\ hypersurface, y/—g-£ = r 2 
sinOdrdOdip. 
Thus 



(3.12c) 



r 



6/7 



and d£ 



\T/ out 

fcim 



1/2 



-2tpXklm 

For t < 0, eq.(3.9) yields the solution 



T- l {l+i)e 



^kim = f- 1/2 \A'J_ 1/2 (r) + B'Y_ 1/2 (f) 



(3.13) 



(3.14a) 



where A' and £?' are integration constants. Here r is defined as 

f = -yJt P X klm {-t + t P a$) (3.146) 

showing that f — > — oo when t — > — oo. 

Using above approximations for Bessel's functions and normalization pre- 
scription, it is obtained that for t — > — oo 



\l/ m 

fcim 



^3 



1/ V 1 (l+2)e +tf . 



(3.15) 



Since $j$m an d both belong to the same Hilbert space, so one write 

k,l,m 
k,l,m 
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As a result, one obtains 



(3.16a) 



where aki m and j3kim are Bogoliubov coefficients satisfying the condition [ 18, 
26-27] 

Wklm? ~ \Pklm\ 2 = 1. (3.16c) 

The in- and out- vacuum states are defined as 



Moreover, 



A)3Jin>=0 = A^|out> 



A out _ jin i/O A*m 

^klm — a klm^klm "r Pklm^klm-> 



(3.17a, 6) 



(3.17c) 



Connecting eqs.(3.11) - (3.17), it is obtained that 



OLklm 

Pklm = ; 

These results yield 
implying 



2 
1 



-i{t+t) 



-i(r+f ) 



/ tj + a 10 tp \ 1/4 / — ix -|- a 10 tp \ 

(: 



*i + a| tpx 1/4 



*i + a? £ P \i/4 



-ti + af tp 



r-( : 



*i + a? *p\ 1/4 



t\ + af tp 
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I 2 - 


-a 10 tp- 


2 


-5tpi 




- *i - 




- ti ^ 



(3.18a, 6) 



(3.19a) 



(3.196) 



2 = 0. 



Eq.(3.19b) shows that when t\ is sufficiently larger than 5tp, \Pkim 
It means that spinless particles will be created in the state < R >= only 
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when t\ < htp. For t\ sufficiently larger than 5tp, there will be no production 
of scalar particles, in this state. 

Anamolous terms emerge when vacuum expectation value of the energy- 
momentum tensor components T^ u is evaluated for $, if dimension of the 
space-time is an even integer. These terms get cancelled for the difference of 
vacuum expectation values T^ u in in- and out-states. Thus T^ u for created 
particles is defined as [28] 

< >=< out|T MI/ |out > — < in|T MI/ |in > . (3.20a) 

for <fi are obtained from the action (2.11) as 

Tfiv = d^d^ - 2 V A[R^ + {(R$*<f>)^ - g^n(M*<S>)} 

-^[V&dpQ - AR 2 $*$] (3.206) 

with its trace 

T = —(1 — 12^^)^$*^$ - 12r?A 2 < R > 3 (3.20c) 
Now trace of for created partcle-antiparticle pairs is obtained as 

< T >=< out|T|out > - < in|T|in > (3.20d) 

Eqs.(3.12b), (3.13), (3.15) and (3.20 d) lead to trace of created particles 
in the state < R >= as 

<Ti(6)>=0. (3.21) 
Case 2: The case of state < R>= \T C 
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Using a(t) from the line-element (2.31) for this state, the equation (3.8) 
is written as 



V 3 /T c - 
2 V 07/ 



kl rn 



a 2 



e -(t-tie)y/T c /&n + jVf 2 



^fcZm = 0, 



(3.22) 



where M 2 = |AT 2 using the definition of M b from eq.(3.2). The equation 
(3.22) yields the solution for t > as 



klm 



7Y 



sinh2na 



1/2 



677 



1/4 _3 

e 4 



(t-t 1Q )^T c /6ri 



X 



where 



and 



1klm& 
2 _ 677 rl 



K*-*io)V T c/ 6r ? 



). 



a = 



-AT 2 - ^ 
4 c 32t/ 



24t/X, 



klm 



a 2 T c 



For £ < 0, the equation (3.22) yields the solution as 



klm 



7]sinh(27ra) 



1/2 



677 



T, 



V4 3,,_ 



e 4 



(t-tl e )y/T c /6 V 



X 



Using the approximation of 

J n (x) 



5(*-*l. 



) v /T c /6 ?7 j 



2 n r(l + n) 

for small x in eqns.(3.23) and (3.24) , for t — > 00 



\l/ out 



r}sinh(2na' 



1/21-677 



(3.23a) 



(3.236) 



(3.23c) 



(3.24) 



1/4 



/ 1 \ rjklm\ 

Vr(l±2ia)A 2 / 



x 
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-(t-ti e )(|Tia) V / ^M 



(3.25a) 



and for t — > — oo 



\T/ m 

klm 



rjsinh2na 



1/2 



6i] 



1/4 / 1 



I 1 \ / 7fc*m \ 



x 



(t-t le )(^Tia)y/T c /6 V 



Using tfj^J, and *jg ro , one obtains 



/3 

i(- =F2ia 



and 



3i 



fcim 



2 

|/3feim| 



(3.256) 



(3.26a) 



(3.266) 



implying that 

(3.26c) 

which ensures creation of spinless particles. Connecting eqs. (3.16c) and 
(3.26b) 



\QZkh 



13 



(3.26d) 



Moreover, the absolute probability for no particle creation is given as 
| < out | in > | 2 = ni a Hm|~ 2 ■ 

klm 

Eqs. (3.26 d) and (3.27) yield 



9 + 16a 2 ' 



(3.27) 



a 2 = 1/4. (3.28) 
Connecting eqs. (3.2), (3.23b) and (3.28), it is obtained that 
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Using the convolution theorem, eq.(3.12b) can be written as 



(3.29) 



®kim(t,r,9,< 



(2n)- l/2 e- ikr 



ro+e 



{y(l - -j)} e^dy]^ klm (t)Yi m (0A) 



(27T) 



-1/2 -ikr_ 



A/7 



(ro + e) 11 /' 



-cos[/c(r + e)\^ k i m {t)Y lm {6, 0). 



(3.30) 



Connecting eqs.(2.27 a), (3.20 d), (3.25) and (3.30) as well as taking aver- 
age over 9 and 0, trace of the energy-momentum tensor for created particles, 
in the state R = (1/2)T C , is obtained at t — t 2e as 



klm 



T m = ^(l/ 27 r)(^) 1/2 S mM(3/4)^(t 2e -t le )] 



X 



(r„ + e) 22 / 7 



cos 2 k(r + e) [(13/16) ^(6?7AT C - 1) - (3/2)^A 2 T ( 



= (2/3n)Big(^) 1/2 S inh[(3/A)^(t 2e - t le )} 



x 



e 8/7 T 



cos 2 A;(r + e)l(13/16)^(6r / AT c - 1) - (3/2)r?A 2 T c 3 
{ro + e) zz ' ' 1 or] 



(3.31) 

Here, summation is done with the help of the Riemann zeta function 
defined as ((s) = D^n -5 ! . Moreover, 



J2cos 2 [k(r + e)] = ^[l-A; 2 (r + e ) 2 + ---] 

klm klm 
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= £i = £(2Z + i) 

klm kl 

= £[4C(-1) + 2C(0)] = -2(4/3)C(0) = 4/3 
k 

as C(— 2m) = 0. Also ((0) = —\ and £(— 1) = — ^ obtained through the 
analytic continuation. 

4. Creation of spin-1/2 particles 

Using invariance of S, given by eq.(2.11) under transformation ip — > + 
Sip, the Dirac equation is obtained as 

(V^ - M f )iJ) = 0, (4.1a) 

where 

M f = aR. (4.16) 

Here, D^ = d^ — T ^ with + 

Y = <l\ (4.2a) 

where (a, // = 0, 1, 2, 3) and are defined through 

= ?U- (4.26) 

Here 77^ are Minkowskian metric tensor components and g^ v are metric tensor 
components in curved space-time. c.c. stands for complex conjugation. 

Dirac matrices 7^ in curved space-time satisfy the anti-commutation rule 

[18] 

{7^, 7"} = 2^ (4.2c) 
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and Dirac matrices 7° in Minkowskian space-time satisfy the anti-commutation 
rule 

{r,i b } = 2 V ab . (4.2d) 



are defined as 



In general, solution of the Dirac equation ip can be written as 



s=±l k 



^ = Y(^Ik,sl° b k,s + ^IIk,sl°d k ,s), 

s=±l k 



where 



and 



With g k l m ,s(f) = f-klm,s{t)- 

In eqs.(4.4) 

ul = (1000), = (0100) 

ul = (0010), = (0001), 



(4.2e) 



(4.3a) 
(4.36) 



ipik,s = [(2^)^(1 - ^) V7 ] lj £fki m ,s(t)yim(0,<l>)e ik - r u a (4.4a) 



fe, s = [(27n7) 1/2 r(l - ^) 1/7 ] _1 E^m,,(t)^ m (^0)e^- r w s . (4.46) 



(4.5a, 6, c, d) 



where the index (T) stands for transpose of the column matrices u± s and 
u± s . In eqs.(4.4), m = —I, • • • , +/; I — 1, 2, 3, • • • and —00 < fc < 00. 
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Using eqs.(3.5) ip is normalized as 



(4.7) 

In eqs.(4.3), b ks (b k!S ) are creation (annihilation) operators for the positive - 
energy particles and d_ ks (d)_ k J are creation (annihilation) operators for the 
negative - energy particles (anti-particles). 

Connecting eqs.(4.2) and (4.4), it is obtained that 

(V/^ - Mf^ip IkiS = 0, (4.8a) 

(i-fDp - M f )rl> IIk>a = 0, (4.86) 
Now using the operator ( — i^D^ — Mfj from the left in eq.(4.8a) as 

( - tYD^ - M f ) (iYD, - M f )^ s = 0, 



it is obtained that [17] 



(□ + (47?)" 1 < R > W < R > 2 )^ IkjS = 0, 



(4.9a) 



where R = r\R and 



□ 



1 d_ 

r^g dx^ 



-99 



d 



dx h 



(4.96) 



Similarly, from eq.(4.7b), it is obtained that 



(□ + (Ar])- 1 <R> +a 2 <R> 2 )i/> IIktS = 0. 



(4.10) 
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When temperature T is not very much below T c one obtains the vacuum 
state < R >= 0, where eqs.(4.7) and (4.8) reduce to 



(4.11a, b) 



Particles are created due to conformal symmetry breaking, which is caused 
by the mass term. In the state < R >= 0, there is no term to break this 
symmetry. So, production of spin-1/2 particles are not possible in this state. 

As given in section 2, one obtains the state < R >= \T C when the 
temperature falls sufficiently below T c . In this state, the universe obeys the 
geometry given by the line-element (2.31) . 

In what follows, creation of spin-1/2 particles is investigated in the model 
given by this line-element. In the background geometry of this model, for 
every k,l,m and s, eq.(4.9) looks like 



Ak.r 



\f ±—f | ({1 | Hjj-_u i 

J klm y s ~r J klm,s '11 2 22 

connecting eqs.(4.4a) and (4.9). Here 



k 2 1(1 + 1) 6r (2r - r ) 



7a 2 r 2 (r — r ) 2 




1 - 



r 



-4/7 



(4.12a) 



M) = {4T])- 1 <R> +a <R> Z 



= ^ + 



(4.126) 

Using the convolution theorem [25] as above and integrating over r, one 
obtains that 

3d • f k 2 r°° 



/ e lk - r dr fan,, + — fkim, s + \~2~ e 



-ik.y 



1 - 



ro 

y 



-4/7 



dy 
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1(1 + 1) 



(2y - r ) 







-4/7 






e 




y- 





6r 
7a 2 i] 



r +ey 2 (y - r ) 



1 - 



-4/7 . 



e- ik -ydy + M 2 f }f klm , s 



0. 



(4.13) 



Details of the evaluation of integrals with respect to y are given in the 
Appendix B. Using these results, in eq.(4.13), one obtains 



J klm,s ~r n \ n J klm,s ' 

2 V 6?7 



X, 



1 lm e -(t-he)y/T c /6v _|_ ]Vf2 



af 



/fcim,* = 0, (4.14a) 



where 



X klm ,s = {l/a 2 rf)e z ' 7 {r + eY' 7 Cos{k{r + e)} 
-yV 2 (r + e)" 1 + !^ V 2 (r + e)" 



2 , + 



(r + e) 2 



Eqs.(4.14) yield the solution for t > 0, 

/«m,. = C' 1 e ± l( t - t -)v^J ±2iQ ( 7WmiS e-^ t - t -)v / ^), 



(4.146) 



(4.15a) 



where Ci is a normalization constant, 



and 



a 



6rj 



f 32rj 



klm,s 



a wT c 



(4.156) 



(4.15c) 



using the definition of M 2 from eq.(4.12b). 
Connecting eqs.(4.5a) and (4.15) 



fe, s = [(2^)^(1 - ^) 1/7 ] _1 X;C7 ie -!(«-«i.)v^x 



28 



^(7«m,.e-5^)V 5 ^)^ m (fl^) e < *'u.. (4.16) 

Normalization of ipik, s is done using the rule (4.7c) at the t = t\ e hyper- 
surface denoted as £ onwards. As a result, C\ is obtained as 



Cl = J±2ia(l 



klm,s i 



(4.17) 



Thus 



1> 



Ik,s 



(2^)^(1 - ^yj x^u^l Vih^ x 



7^se^ (t - t - ) ^M)^ m (^ 0) e ^ Ms . (4.18) 



Similarly 



IIk,S 



(2^)^(1 - ^) i/ v i Ek ± ^(7-^)r i ^ |( " ie) ^x 



J ±2M (7-^, s e^(^-)v / ^)^ m (^0)e- fe -^. (4.19) 
Asymptotics of these solutions, when t — > oo, yield 



out 
Ik,s 



(2^/V(l - ^) 1/ V 1 E|^(7H m , s )r 1 e-t^)v / ^x 



r(l-2ia)] 2l V i ^-^)v / ^F Zm (^0)e^ Ms . (4.20a) 



and 



7 /,out 

riik,s 



" 1 e -!(*-*ie)V Tc / 6? ?x 



HI -2ia 



(2 mj ) I / 2 r(l-^) / ']'E|j J3ia ( 7 -H mi( ) 

1 (2z|l^)-^ e ^«(^i«)v^^ ro (0,0) e -^ 8 . (4.206) 
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For t — > — oo, one obtains 

= [(2^)^(1 - ^) 1/7 ]"Ek^(7^)r 1 e^ (t - tle) v / ^x 



r(l + 2m)] 1 2 V^-^v^7^ m (0, 0)e ifc - r Ms . (4.21a) 



and 



= [(2^)^(1 - ^) 1/ V 1 E|^(7-H m , s )p 1 ef^)v^x 



r(l + 2ia)] '(Iz^) 2ia e--(^e)V^Y ;m (^0)e- ife -^ s . (4.216) 

In eqs.(4.20) - (4.21), T(x) stands for the gamma function of x. 
In the in- as well as out- region, 
the decomposed form of ip can be written as [4] 



j Ti(-k,s) 



s=±l fc 

= E E(^M^°(M + ^-C-sV'°/(-fc,-s)) ' 



s=±l fc 

(4.22a, 6) 

as both in- and out-spinors belong to the same Hilbert space. The in- and 
out- vacuum states are defined as 



&£J in >= cT M | in >= 



(4.23a, 6) 



and 



6^ | out >=(T£ s |out >=0 



(4.23c, d) 
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Bogoliubov transformations are given as [3,4] 



tout 
°k,s 

i fout 
°k,s 

dlt-s 
d°-l-s 



k,—s 



k^s k^s > k 

b™ s a kyS + d ] ll_ s f3 Ks 

k,s k,s ik.s k. s 



Connecting eqs.(2.22)-(2.24), one obtains [22] 

|afcr + |/5fcr = EK-| 2 + IW 2 = 1 ' 



(4.24a, b, c, d) 



(4.25) 



where 



and 



Connecting eqs. (4.20a), (4.21a) and (4.26a) 



a k . s = (l/2)[r(l-^)]^(^)" 4ja 



(4.26a) 



(4.266) 



J±2ia{lklm,s) 



sinh(2na) 
(Ana) 



r(l-2ia) 



(4.27a) 



Similarly eqs.(4.20b), (4.21b) and (4.26b) 

P k . s = (l/2)[T(l + 2ta)Y\^^y' ia \j ±2ia (lki my s) 



sinh(2ira) 
(Ana) 



T(l + 2ia) 
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Pkjak.8 2 = 1- (4-29) 



(4.276) 

Conditions (4.25) and eqs.(4.27) yield an useful result 

Ksl 2 = \Pk.s\ 2 = \ (4.28) 

suggesting that a should be very small, using the symmetry k — > —k caused 
by symmetry of the model under r — ■> — r. 

The relative probability of creation of a particle - antiparticle pair is given 

as 

^klm,s 

Absolute probability of the creation of particle - antiparticle pairs requires 
the total probability of creating 0, 1, 2, • • • pairs to be unity [5], which means 
that 

N k lm,s(l + U klm y s + ^klm,s H ) = ^ 

where Nki m ,s is the probability of creation of no pair of particle-antiparticle. 
Using eq.(4.29), it is obtained that 

Nhlm >i = = 0. 

It shows that probability of vacuum to remain vacuum is 

| < out|in > | 2 = II N kim,s = 0. (4.30) 

klm,s 

implying certainity of creation of particle-antiparticle pairs. 

The action S, given by eq.(2.11), yields components of energy - momen- 
tum tensor for ip field as [22,27] 

1 5_S 
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T 

-ori(D^D v - g^n^tp + c.c. 



(4.31) 



Connecting eqs. (3.20a) and (4.31), energy density of created particles in 
the state R = 1/2T C is obtained as 

p =< T ° > = < in| - i$7°d - a^R^ip + 2oTh\)T\) + c.c. |in > 

— < out | — i-^7°(9o — o-r/RQipip + 2aR^ + c.c. | out > 
= < in | — 2^7°c\) + l/2aT c ipi/j + c.c. |in > 

- < out | - m/>7°<9 + l/2aT c ^/> + c.c. | out > 

(4.32) 

using = 3d /a = T c /2r]. 

In- and out- ip, given by eqs. (4.20) and (4.21), recast eq.(4.32) as 

g a 8/7 

P = E —7 ri2277 cos2 [ fc ( r o + e)]sm/i[(3/2)(t - t u )jT c /6rj\ (4.33) 

using convolution theorem for ip as given in the preceding section. 

Similarly trace of stress - tensor for created spin-1/2 particles is given as 

< T >=< in|T^|in > - < out|T M „|out >= 0, (4.34) 

using In- and out- solutions for ip, given by eqs. (4.20) and (4.21). 
Eq.(4.28) shows number density of created spin-1/2 particles as 
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IAI 2 = \ (4.35) 
yielding energy density of created particles also as 

P = \M S . (4.36) 

|/9fc| 2 gives number of created particles per unit volume during the time period 
(^20 — ^io) — 1-38 x 10 7 tp, so evaluating p, from eq.(4.32) and comparing with 
the same from eq.(4.36), it is obtained that 

8/7 

M, - 54.87 x 10" (rQ + e)22/7 . (4.37) 

Here also summation is done with the help of the Riemann zeta function 
defined in the earlier section. 

4. Conclusions 

Thus, in this part of the series, production of spinless and spin- 1/2 par- 
ticles are discussed in the first two phases of the universe. Contribution of 
these particles to later development of the universe to third and fourth phases 
of the universe will be discussed in paper II. 
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Riccion and Graviton 



From the action 



S = Jdxd y^-g { 4 +D) [ — 



+ a (A+D)R(4+D) + 



7(4+D)(-R(4+D) jy^2y n (D+i)R 2 (D+A))_ > ( A -l) 



the gravitational equations are obtained as 



M( 2+D ) 1 (1) (2) 

— — {RmN — -9MnR(A+D)) + &(A+D)H MN + 7(4+d)^n = 0, (A2a) 

where 



-^mtv — 2R.MN — ^9mn^(a+d)R(a+d) — 29mnR\a+d) + 2R(4+d)Rmn, {A. 2b) 
and 

HmN — 3R^MN - ^9MN U (A+D)R 2 (A+D) ^-^y{~-^9MN D (A+D)R 2 4+D) 

+2D( 4+£ ))i?( £ ) + 4)i? M jv + R 2 ,mn) ~ 29mnR\a+d) + 3R 2 (a+d)Rmn- (A 2c) 

Taking ^Miv = Vmn+Iimn with ttmjv being (4+D)-dimensional Minkowskian 
metric tensor components and Hmn as small fluctuations, the equation for 
graviton are obtained as 



u (A+D)h M N = (A3) 

neglecting higher-orders of h. 
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On compactification of M 4 <g> S D to M 4 , eq.(A.3) reduces to the equation 
for 4-dimensional graviton as 

for the space time 



dS 2 = g^dx^dx v - p 2 [d6l + sin 2 6id6 2 2 ^ \-sin 2 9 1 ■ ■ ■ sin 2 Q (D ^ x) dQ 2 D . (A5) 

The 4-dimensional graviton equation (A. 4) is like usual 4-dimensional 
graviton equation ( the equation derived from 4-dimensional action) only for 
I = 0. Thus the massless graviton is obtained for I = only. 

As explained, in section 2, the trace of equations (A. 2) leads to the riccion 
equation 

[□ + l -^R + m\ + ±R 2 ]R + = 0, (A6a) 

where 



£ = 



A 



D 



2(D + 3) 
(D + 2)\V D 

1 



rj 2 XR 



D 



+ 



2(D + 3) 2 1 D 



A(D + 3)a 

(D + 2)\M^ 2+D ^V D 



-v 



16tt 



+ 



DR 2 D 



4(0 + 3) + 
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The graviton has 5 degrees of freedom (2 spin-2 graviton, 2 spin-1 
gravi- vector (gravi-photon) and 1 scalar). The scalar mode / satisfies the 
equation 

P 

Comparison of eqs.(A.6) and (A. 7) show many differences between scalar 
mode / of graviton and the riccion (R) e.g. £, A and d , given by eqs.(A.6b,c,d,e), 
are vanishing for /, but non- vanishing for R. Eq.(A.7) shows (mass) 2 for / 
as 

= ifl + s- i) 

whereas (mass) 2 for R, given by eq.(A.6c), depends on G^ + d,Vd and Rd 
(given in section 2). m 2 = for / = 0, but m 2 ^ can vanish only when gravity 
is probed upto ~ lCT 33 cm. As mentioned above, so far, gravity is probed 
only upto 1cm. Thus (mass) 2 of riccion does not vanish. 

So, even though, / and R are scalars arising from gravity, both are dif- 
ferent. Riccion can not arise without higher-derivative curvature terms in 
the gravitational action, but graviton can be obtained even from Einstein- 
Hilbert action. 

1 d f , „„ d x „„ d 



where X M are locally inertial co-ordinates and n^ are Minkowskian metric 
components. It shows that the scalar like operator □ has the same role on 
R as it is for other scalar fields due to principle of equivalence. 
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Appendix B 



In eq.(3.6), using the convolution theorem the function e tk ' r 
can be written as 



I _ ra 



-4/7 







—4/7 f°° 

= r] -i e Lk.r dye ~ik.y 




-4/7 




L r _ 


Jr +e 







where 77 has dimension of length, as used in section 2, for dimensional correc- 
tion. Here in the integral dy introduces a quantity of length dimension which 
is compensated by 77. 

Evaluation of the integral in eq.(B.l), is done as follows. 



/ dy 

Jro+e 



-ik.y 



\ — IS. 

y 



4/7 



= E 



r 0+ e y (y-r )^/^ =0 

(-ik) n 1 



^ n\ T(u + 4/7)J n)+ e 

x u,-3/7 e -x(y-r ) dx 

-ik) n 1 



dyy n+4/7 (y - r ) u x 



E 



^ n! I> + 4/7)70 
[J" y n+Ap {y-r,Te- x ydy\ 



dxx"- 3/7 e xro) x 



- E 



-ik) 1 



1 



„ „ n\ r(^ + 4/7)7o 
1 



(r + e 



( (n+4/7) e o; 



+ — + 

x ear 



dxx"- zr ' 'e~ xe x 
uj uj{uj — 1) 



e 2 x 3 



+ 



(B.2) 



neglecting terms containing (r + e) 1 compared to terms containing e 1 . 
Performing further integration in eq.(B.2), one obtains that 



poo 

/ dy 

Jro+e 



-ik.y 



y 



4/7 



- E 



' n! I> + 4/7) 



(r + e)( n+4 / 7 ) x 
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,3/7 



T(uj - 3/7) + uT(lu - 10/7) + 



u(u - l)V(u - 17/7) + 



E- 



-i£;) n 1 



(r + e 



,(n+4/7) 3/7 



= (r + e) 4/7 e 3/7 cosk(r + e), 



(5.3) 



when uj — > oo. Here, symmetry under fc — > —A; caused by symmetry of the 
model under r — > — r is used. 
The same procedure yields 

/°° dye- ik *y- 2 [l - -1 ^ = (r + e)- 10 / 7 e 3 / 7 cosA;(ro + e), (5.4) 



/ dye- iky y- 3 / 7 [y - r ]- 18/r = (r + e)- 3 /V n / 7 cosA;(r + e), (5.5) 

/"OO 

/ dye lk - y y- 4/7 [y - r ] 4/7 = (r + e)- 4/7 e n / 7 cosA;(ro + e), (5.6) 

/ dye- ik - y y- 10 ^[y - r ]~ 18/7 = (r + e)~ w / 7 e' 11 ' 7 cosk{r + e). (5.7) 
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